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Abstract. We present XMM-Newton MOS imaging and PN timing data of the faint millisecond pulsars PSR J075 1 + 1807 and 
PSR J1012+5307. We find 46 sources in the MOS field of view of PSR J0751 + 1807 searching down to an unabsorbed flux 
limit of 3 xlO"" ergs cm"' s"' (0.2-10.0 keV). We present, for the first time, the X-ray spectra of these two faint millisec- 
ond pulsars. We find that a power law model best fits the spectrum of PSR J075 1 + 1807, r=1.59±0.20, with an unabsorbed 
flux of 4. 4x10"''' ergs cm"-s"' (0.2-10.0 keV). A power law is also a good description of the spectrum of PSR J1012+5307, 
r=1.78±0.36, with an unabsorbed flux of 1.2xl0"'^ergs cm"^s"' (0.2-10.0 keV). However, a blackbody model can not be ex- 
cluded as the best fit to this data. We present some evidence to suggest that both of these millisecond pulsars show pulsations 
in this X-ray band. We find some evidence for a single broad X-ray pulse for PSR J075 1 + 1807 and we discuss the possibility 
that there are two pulses per spin period for PSR J1012+5307. 
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1. Introduction 

X-ray emission from millisecond pulsars (MSPs) is thought 
to be from: charged relativistic particles accelerated in the 
pulsar magnetosphere (non-thermal emission indicated by a 
hard power- law spectrum and sharp pulsations); and/or thermal 
emission from hot polar caps; and/or emission from a pulsar 
driven synchrotron nebula; or interaction of relativistic pulsar 
winds with either a wind from a close companion star or the 
companion star itself (see e.g. Becker & Aschenbach 2002, 
for a more thorough review). Until now. X-ray pulsations and 
spectra have often not been observable for faint MSPs e.g. 
PSR J0751 + 1807 (Becker et al. 1996) and PSR J1012+5307 
(Halpern & Wang 1997). Thus it has been difficult to discrimi- 
nate between competing neutron star models. However, taking 
advantage of the large collecting area of XMM-Newton (Jansen 
et al. 2001), it is becoming possible to observe not only the 
X-ray spectra but also put a limit on the presence of X-ray pul- 
sations of these faint MSPs. 

We have observed two faint MSPs with XMM-Newton. PSR 
J075 1 + 1807 was first detected in an EGRET source error box, 
in September 1993 (Lundgren et al. 1993), using the radio 
telescope at Arecibo. Lundgren et al. (1995) combined the 
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mass function, eccentricity, orbital size and age of the pulsar, 
determined from radio data, to predict the expected type of 
companion star to the millisecond pulsar. They proposed that 
the secondary star is a helium white dwarf, with a mass be- 
tween 0.12-0.6Mo, in a 6.3 hour orbit with the pulsar. They 
determined a 3.49 ms pulse period, but from the period deriva- 
tives, the spin down energy indicated that the pulsar is not the 
source of the y-rays that were originally detected by EGRET. 
PSR J075 1 + 1807 was subsequently detected in the soft X- 
ray domain by Becker et al. (1996), using the ROSAT PSPC. 
However, there were too few counts detected to fit a spectrum 
or detect pulsations. Using the HI survey of Stark et al. (1992), 
they deduced an interstellar absorption of 4x 10^*' cm"^. Then 
using their estimated counts and assuming a powerlaw spec- 
trum of dN/dE oc E^^-^, they determined an unabsorbed flux of 
1 X 10"'^ erg cm"^ s"' (0. 1-2.4 keV). The corresponding X-ray 
luminosity was then calculated to be = 4.7 x 10-" erg s"', 
for a distance of 2 kpc. The distance they used was calculated 
using the radio dispersion measure and the model of Taylor & 
Cordes (1993) for the galactic distribution of free electrons. 

PSR J1012+5307, a 5.26 msec pulsar, was discovered by 
Nicastro et al. (1995) during a survey for short period pulsars 
conducted with the 76-m Lovell radio telescope. From the dis- 
persion measure they found a distance of 0.52 kpc. Callanan et 
al. ( 1 998), and references therein, confirmed a white dwarf sec- 
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ondary of mass O.I6+O.O2M0 in a 14.5 hour circular orbit with 
the pulsar. Halpem (1996) associated a faint (L^: «2.5x 10^" 
ergs s~', 0.1-2.4 keV) X-ray source detected with the ROSAT 
PSPC (80+24 photons), with the radio MSP PSR J1012-H5307. 
However, the number of photons was insufficient to determine 
a spectrum or any pulsations. 

In this work we present the X-ray spectrum of both PSR 
J0751-I-1807 and PSR J1012-I-5307 for the first time. We also 
present some evidence for X-ray pulsation from both of these 
faint pulsars and we compare their nature with other millisec- 
ond pulsars, e.g. PSR J0437-4715 (Becker & Triimper 1993; 
ZavUn& Pavlov 1998; Zavlin et al. 2002). 

2. Observations and data reduction 

PSR J0751-H1807 was observed by XMM-Newton on 2000 
October 1. The observations spanned 38 ksecs (MOS cam- 
eras) and 36.8 ksecs (PN camera), but a solar flare afifected ap- 
proximately 8 ksecs of these observations. The second of our 
faint MSPs, PSR J1012-I-5307, was observed on 2001 April 19. 
The MOS observations lasted 20.8 ksecs and the PN observa- 
tions 19.2 ksecs. However the whole observation was strongly 
aff"ected by a solar flare. The MOS data were reduced us- 
ing Version 5.4.1 of the XMM-Newton SAS (Science Analysis 
Software). However, for the PN data we took advantage of 
the development track version of the SAS. Improvements have 
been made to the oal (ODF (Observation Data File) access 
layer) task (version 3.106) to correct for spurious and wrong 
values, premature increments, random jumps and blocks of 
frames stemming from diff'erent quadrants in the timing data 
in the PN auxiliary file, as well as correcting properly for the 
onboard delays (Kirsch et al. 2003a). Indeed several spurious 
jumps were corrected using this version. We have verified that 
this version improved our liming solution using the pulsar PSR 
J0218+4232 (see Webb etal. 2003). 

We employed the MOS cameras in the full frame mode, 
using a thin filter (see Turner et al. 2001). The MOS data 
were reduced using 'emchain' with 'embadpixfind' to detect 
the bad pixels. The event lists were filtered, so that 0-12 of 
the predefined patterns (single, double, triple, and quadruple 
pixel events) were retained and the high background periods 
were identified by defining a count rate threshold above the low 
background rate and the periods of higher background counts 
were then flagged in the event list. We also filtered in energy. 
We used the energy range 0.2-10.0 keV, as recommended in 
the document 'EPIC Status of Calibration and Data Analysis' 
(Kirsch et al. 2002). The event lists from the two MOS cameras 
were merged, to increase the signal-to-noise. 

The PN camera was also used with a thin filter, but in tim- 
ing mode which has a timing resolution of 30yus (Striider et 
al. 2001). The PN data were reduced using the 'epchain' of 
the SAS. Again the event lists were filtered, so that 0-4 of the 
predefined patterns (single and double events) were retained, 
as these have the best energy calibration. We also filtered in 
energy. The document 'EPIC Status of Calibration and Data 
Analysis' (Kirsch et al. 2002) recommends use of PN tim- 
ing data above 0.5 keV, to avoid increased noise. We used the 
data between 0.6-10.0 keV as this had the best signal-to-noise. 



The times of the events were then converted from times ex- 
pressed in the local satellite frame to Barycentric Dynamical 
Time, using the task 'barycen' and the coordinates derived 
from observing the pulsar 4-8 times per month, over years, with 
the radio telescope at Nan9ay, France, by one of us (IC), for 
PSR J0751 + 1807. For PSR J1012+5307 we used the results 
from the high-precision timing observations taken since 1996 
October, using the 100-m Efl'elsberg radio telescope and the 
76-m Lovell telescope by two of us (MK and OL). 

3. PSR J0751+1807 

3.1. Spectral analysis 

We extracted the MOS spectra of PSR J075 1 + 1807 using an 
extraction radius of ~l' and rebinned the data into 15 eV 
bins. We used a similar neighbouring surface, free from X-ray 
sources to extract a background file. We used the SAS tasks 
'rmfgen' and 'arfgen' to generate a 'redistribution matrix file' 
and an 'ancillary response file' . We binned up the data to con- 
tain at least 20 counts/bin. The PN data were extracted in a sim- 
ilar way, following the usual XMM-Newton timing data proce- 
dure. The data in the RAWY direction were biimed into a single 
bin and the spectrum was extracted using a rectangle of 1 x 3 
pixels, which included all the photons from the pulsar, see e.g. 
Kuster et al. (2002). The background spectrum was extracted 
from a similar neighbouring surface, free from X-ray sources. 
Again we used the SAS tasks 'rmfgen' and 'arfgen' to generate 
a 'redistribution matrix file' and an 'ancillary response file' . We 
then used Xspec (Version 11.1.0) to fit the spectrum. We tried 
to fit simple models to the combined PN and MOS spectra. We 
find the model fits as given in Table 1 for the spectrum between 
0.2-10.0 keV, when the Nh was frozen at 4 xlO^^cm"^, (see 
Becker et al. 1996). We find that the best fitting spectrum is a 
single power law, with a similar photon index to the X-ray spec- 
trum of other millisecond pulsars e.g. PSR B1821-24 (Saito et 
al. 1997). The spectrum of PSR J0751 + 1807, plotted wifli the 
power law fit, can be seen in Fig. 1. Allowing the Nh to vary, 
gives values compatible with the above values. We determine 
an unabsorbed flux of 4.4 xlO-^'' ergs cm'^s"! (0.2-10.0 keV). 

3.2. Timing analysis 

We have reduced and analysed the timing data in the same way 
as the timing analysis that we carried out on the MSP PSR 
J021 8+4232 (Webb et al. 2003), which is a MSP that is known 
to show pulsations in both the radio and X-rays. We corrected 
the timing data for the orbital movement of the pulsar and the 
data were folded on the radio ephemeris, see Table 2, taking 
into account the time-delays due to the orbital motion. We used 
the data between 0.6-7.0 keV, as we found that the majority of 
the emission from PSR J075 1+1807 was in this energy band 
(see Sect. 3.1) and thus the signal-to-noise in this band was 
the best. We tested the hypothesis that there is no pulsation in 
the MSP PSR J075 1+1807. We searched frequencies at and 
around the expected frequency and we found the largest peak 
in the xl versus change in frequency from the expected fre- 
quency at the expected value, see Fig 2 where we have taken 
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Table 1. The best model fits to the MOS spectra of the source 52, and the MOS and PN spectra of PSR J0751 + 1807 and PSR J1012+5307 . 
For PSR J0751+1807 the interstellar absorption ( N/,) was fixed at 4 x 10^" cm ^ For PSR J1012+5307 the was fixed at 7 x 10" cm ^ 
For source 52, z was fixed at 0.255. All the errors shown are Icr confidence limits. The fluxes and luminosities are given for the 0.2-10.0 keV 
band. 



Object 


Spectral 


Nh (cm-^) 


kT 


Photon 


Abundance 


xl 


dof 


Flux(xlO-'^') 


Luminosity 




model 


(X 10^2) 


(keV) 


Index 








(ergs cm"^s"') 


(ergs s"') 


PSRJ0751+1807 


Power law 


0.04 




1.59+0.20 




1.33 


14 


0.44 


2.1xlff" 




Bremsstrahlung 


0.04 


12.36+12.31 






1.52 


14 








Blackbody 


0.04 


0.32+0.04 






2.21 


14 






Source 52 


Power law 


0.13+0.06 




1.46+0.13 




0.86 


32 


2.0 


1.9xl0« 




Bremsstrahlung 


0.09+0.05 


17.31+10.42 






0.87 


32 








Raymond Smith 


0.07+0.04 


17.63+12.78 




0.16+1.20 


0.89 


31 






PSRJ1012+5307 


Power law 


0.007 




1.78+0.36 




1.27 


9 


1.2 


3.9x10^" 




Bremsstrahlung 


0.007 


1.95+1.45 






1.29 


9 








Blackbody 


0.007 


0.26+0.04 






1.37 


9 







Table 2. Ephemeris of PSR J075 1 + 1 807 from the Nangay radio timing 
data. 




Energy (keV) 

Fig.l. The combined MOS and PN spectrum of PSR J0751 + 1807 
fitted with a power law model. The fit parameters can be found in 
Table 1. 



Parameter 


Value 


Right Ascension (J2000) 


07" 51"" 09^56312 


Dechnation (J2000) 


18° 07' 38"590620 


Period (P) 


0.003478770781560571 s 


Period derivative (P) 


0.726912 xlO-^o s s"' 


Second period derivative (P) 


2.48928 xlO-^" ss-2 


Frequency (v) 


287.457858811 Hz 


Frequency derivative (y) 


-6.0066207 xlO-'«Hzs-' 


Second frequency derivative (v) 


-2.0569423 xlO'^^Hzs-^ 


Epoch of the period (MJD) 


49301.5 


Orbital period 


22735.664643 s 


a.sin i 


0.39660728 


Eccentricity 


0.0000566981 


Time of ascending node (MJD) 


49460.430540 



the resolution of the data (n) and plotted 4n bins in the range 
-l.TSxlO-"* < Af < 1.75 x 10 1 Testing the significance of 
the peak (Buccheri et al. 1985), we find that it is significant 
at 1.7cr. However, as this is the largest peak when searching 
700 frequencies about the expected frequency and it falls at 
the expected value, we tried to fold the data on this frequency. 
The folded lightcurve (0.6-7.0 keV), counts versus phase, is 
shown in Fig 3. We find one broad pulse per period. Fitting 
the lightcurve with a Lorentzian (as Kuiper et al. 2002a) we 
find that the FWHM of the pulse is =0.31 1+0.1, centred at 
phase 0i=O.38±O.O4 (errors are 90% confidence). Fitting with 
a Gaussian gives similar results. Using a test (Buccheri et 
al. 1983), which is independent of binning, we determine a 
value of 5.5, which corresponds to a probabihty that the pulse- 
phase distribution deviates from a statistically flat distribution 
of 0.94. The pulsed percentage is 52±8%. 

Taking the 261 ±60 counts above the background in our 
observation, we are also able to analyse the variations in the 
lightcurves from different energy bands. We chose the energy 
bands: 0.6-1.0 keV; 1.0-2.0 keV; and 2.0-7.0 keV. The data 
were folded as before and biimed into 5 bins, where two cy- 



cles in phase, are presented for clarity. These lightcurves can be 
seen in Fig. 4. As can also be seen from the spectral fitting, see 
Sect. 3.1, PSR J075 1 -I- 1 807 emits mostly at lower energies. The 
lightcurves do not vary dramatically, although the peak may be 
slightly narrower at lower energies. The pulsed percentage is 
possibly less at lower energies. We find 32±16% for each of 
the energy bands: 0.6-1.0 keV; and 1.0-2.0 keV and 55±23% 
for the 2.0-7.0 keV band. 

3.3. The MOS field of view 

We have detected 46 sources in total using the task 'emlde- 
tect', with a maximum detection likelihood threshold of 4.5cr 
and ignoring those sources not found on both cameras (un- 
less they lay outside the FOV or were in a chip gap). These 
sources can be seen in Fig. 5 and details such as their position, 
count rates and likehhood of detection can be found in Table 3. 
PSR J0751-I-1807 is labelled as source 47. RO SAT detected 10 
sources in the XMM-Newton FOV (Becker et al. 1996). We 
calculate an unabsorbed flux hmit of our observation of 3.0 
xlO-'^ergs cm-^s-i (0.2-10.0 keV), using PIMMS (Mission 
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A f 



Fig. 2. xt versus change in frequency from the expected pulsation 
frequency (shown as the solid vertical line at A f = 0.0) for PSR 
J075 1+1807. 



1.0 

Phase 



Fig. 3. Lightcurve of PSR J075 1 + 1807 folded on the radio ephemeris 
and binned into 6 bins, each of 0.58 msecs. Two cycles are shown for 
clarity. A typical ±lcr error bar is shown. The dashed line shows the 
background level, where the error bar represents the ilcr error. 



Count Rate Simulator) Version 3.3a, with a power law, T = 2 
(as Hasinger et al. 2001). 

The brightest source in the field of view is source 52, 
which was not detected by ROSAT, where the deepest ob- 
servation had a detection limit of 1.5 XlO~^* ergs cm~^s~^ 
(0.2-10.0 keV). We detect this source with an unab- 
sorbed flux of 2.0 xlO''^ ergscm-^s"^ (0.2-10.0 keV), 
more than 13 times brighter than the ROSAT detection 
Umit. This source has been observed as a part of the 
AXIS (An XMM-Newton International Survey) project (see 
http://www.ifca.unican.es/~xray/AXIS/and also Barcons et al. 
2002). The source identified as the optical counterpart of our 
source 52 in the i-band photometry is elongated and the spec- 
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Fig. 4. Lightcurves of PSR J075 1 + 1 807 from different energy bands, 
folded on the radio ephemeris and binned into 5 bins. Two cycles are 
shown for clarity, (a) 0.6-1.0 keV; (b) 1.0-2.0 keV; (c) 2.0-7.0 keV. A 
typical ±lcr error bar is shown. 
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Fig. 5. The MOS field of view. The sources detected with a likeli- 
hood >4.5cr are encircled and labelled with a number. These num- 
bers correspond to the source identification numbers in Table 3. PSR 
J0751+1807 is labelled 47. 



trum shows calcium H and K lines and the G-band. From the 
spectral line shifts, a redshift of z=0.255 has been determined 
by the AXIS group and they propose that this source could be 
a galaxy. We extracted the MOS spectra of source 52 in the 
same way as described in Section 3.1. We have fitted the spec- 
trum with the extra-galactic origin in mind and the results can 
be found in Table 1. The data fitted with a power law model 
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Table 3. X-ray sources in the PSR J0751+1807 MOS field of view. 
Information includes the source identification number, right ascension 
and declination of the source, count rate and detection likelihood. The 
detection likelihood value is the value given by 'emldetect' for de- 
tection in the two MOS cameras, corrected for the coding error in 
versions 4.1 1.15 and earlier. 



Src 


R.A. (2000) 


dec. (2000) 


Counts/s 


Like- 


ID 


h m s 




xlO-3 


lihood 


3 


7 50 54.1 


17 57 59.3 


5.6±0.8 


62 


4 


7 50 54.8 


17 58 27.9 


5.6±0.8 


63 


5 


7 51 12.9 


17 58 44.3 


4.4±0.6 


53 


13 


7 51 9.7 


18 1.04 


3.3±0.5 


40 


17 


7 51 46.8 


18 36.0 


6.4+0.8 


99 


19 


7 51 44.0 


18 1 7.05 


3.8±0.7 


32 


20 


7 51 13.8 


17 57 21.3 


3.7+0.7 


23 


23 


7 50 45.1 


18 1 30.8 


3.5±0.6 


31 


24 


7 50 25.2 


18 1 48.4 


6.6±1.1 


40 


25 


7 51 26.9 


18 2 0.08 


6.4±0.6 


134 


26 


7 50 20.9 


18 2 25.6 


5.3±1.0 


34 


27 


7 50 57.4 


18 3 19.2 


2.8+0.4 


39 


28 


7 51 1.6 


18 3 20.7 


3.7±0.5 


65 


29 


7 51 35.6 


18 3 29.6 


4.9±0.6 


64 


31 


7 50 36.6 


18 3 28.8 


3.0±0.7 


18 


33 


7 51 7.9 


18 4 21.4 


3.6±0.4 


67 


34 


7 50 37.9 


18 5 2.7 


3.3+0.6 


28 


36 


7 50 52.3 


18 5 19.1 


2.6±0.4 


32 


37 


7 50 38.0 


18 5 40.8 


4.0±0.6 


48 


38 


7 50 48.1 


18 5 55.1 


3.3+0.5 


60 


39 


7 51 40.4 


18 6 12.9 


9.6±1.1 


106 


40 


7 50 49.2 


18 6 39.9 


2.8±0.4 


43 


41 


7 52 1.1 


18 6 42.6 


5.3+0.9 


38 


42 


7 51 48.5 


18 6 44.6 


8.3+0.8 


141 


43 


7 50 52.7 


18 6 55.7 


6.9+0.6 


202 


44 


7 51 13.6 


17 56 50.7 


6.1+0.9 


56 


45 


7 51 37.2 


18 7 20.4 


9.6+0.7 


255 


47 


7 519.1 


18 7 36.3 


8.9+0.6 


327 


51 


7 51 4.2 


18 8 45.0 


12.6+0.7 


466 


52 


7 51 17.8 


18 8 56.4 


36.3+1.2 


2017 


55 


7 51 51.5 


18 9 59.2 


4.2+0.6 


44 


56 


7 50 46.5 


18 9 45.2 


4.6+0.7 


56 


59 


7 50 41.7 


18 10 26.7 


3.0+0.5 


37 


60 


7 51 27.3 


18 10 40.2 


3.0+0.5 


41 


62 


7 51 14.9 


18 11 30.0 


3.1+0.4 


56 


63 


7 50 40.1 


18 12 14.3 


17.2+1.2 


355 


64 


7 51 17.5 


18 12 2.9 


2.7+0.4 


47 


65 


7 51 56.1 


18 13 55.7 


5.4+1.0 


37 


66 


7 51 54.1 


18 13 54.9 


8.6+1.2 


73 


67 


7 51 20.3 


18 14 0.7 


6.6+0.6 


117 


68 


7 51 37.5 


17 55 46.9 


6.9+1.1 


53 


71 


7 51 13.5 


18 15 33.4 


3.7+0.6 


40 


72 


7 50 43.2 


18 15 41.3 


4.4+0.8 


28 


73 


7 51 36.8 


18 16 0.5 


6.3+0.8 


96 


75 


7 51 24.2 


18 17 43.5 


5.2+0.8 


38 


76 


7 51 17.5 


17 54 10.7 


6.1+1.0 


41 



can be found in Fig. 6. We investigated both models with the 
redshift frozen at 0.255 and also models with the redshift as 
a free parameter. We found that the fits for the frozen and the 
variable redshift (z) were the same within the error bars, hence 
we present only one set of results, those where the redshift 




0.5 1 2 5 

Energy (keV) 



Fig. 6. The spectrum of the source 52, in the PSR J075 1-1-1807 field of 
view, fitted with a power law model. The fit parameters can be found 
in Table 1. 

was fixed at 0.255. We find a reasonably hard X-ray spectrum, 
which can not be fitted with a black body. We have also inves- 
tigated a temporal variation, binning the data into Iks and 5ks 
bins. However, we find no significant variation during our ob- 
servation. From the X-ray spectral fitting and short time scale 
temporal analysis, we can not support nor disregard the galaxy 
hypothesis. 

4. PSR J1 01 2+5307 

4.1. Spectral analysis 

Due to the observation of PSR J1012+5307 being affected by a 
solar flare during the whole of the observation, we detect only 
one source in the whole of the field of view. This source is at the 
centre of the field of view, with the same coordinates as those 
of PSR J1012-I-5307 and is therefore supposed to be the MSP 
PSR J1012-I-5307. We have not carried out any filtering for pe- 
riods of higher background, as the whole observation would be 
included. Thus the signal to noise is poor for this observation. 

We extracted the spectra in the same way as described in 
Sect. 3. We tried to fit simple models to the combined PN and 
MOS spectra. We find the model fits as given in Table 1 for the 
spectrum between 0.2-10.0 keV, when the Nh was frozen at 7 
xlO'^cm"^ (see Snowden et al. 1994). We can not discriminate 
which of these fits is the best, however we present the spectrum 
of PSR J1012-I-5307 witii a power law fit in Fig. 7. Allowing 
the Nh to vary, gives values compatible with the above values. 
We determine an unabsorbed flux of 1.2 xlO"'^ergs cm"^s~' 
(0.2-10.0 keV). 

4.2. Timing analysis 

We corrected the timing data, as before, for the orbital move- 
ment of the pulsar and the data were folded on new radio 
ephemerides calculated to be correct for our observation (see 
Table 4), taking into account the time-delays due to the orbital 
motion. We used the data between 0.6-5.0 keV, as the signal-to- 
noise was best in this band. Also from the spectral fitting (see 
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Fig. 7. The combined MOS and PN spectrum of PSR J1012+5307 
fitted with a power law model. The fit parameters can be found in 
Table 1. 



Sect. 4.1) we found that the majority of the emission from PSR 
J1012+5307 was in this energy band. We tested the hypoth- 
esis that there is no pulsation in the MSP PSR J1012+5307. 
The largest peak in the xl versus change in frequency from the 
expected frequency is at Av ^ 8 x lO^^'s"', see Fig 8. This 
implies a Af/f (=AP/P) of ~4 xlO"**, similar to the value that 
we found for PSR J0218+4232 (Webb et al. 2003) (~1 xlO"**) 
and approaching the value of < 10 determined from the anal- 
ysis of two revolutions of data of the MSP the Crab (Kirsch et 
al. 2003 a,b). It is also well inside the resolution of this dataset 
{~l/T„i,s), which is 7 xlO^^s"^ thus we can conclude that the 
data reduction and analysis made to the dataset are reliable. 
Testing the significance of the peak (Buccheri et al. 1985), 
we find that it is significant at 3cr. We folded the data on the 
radio period given in Table 4. The folded lightcurve (0.6-5.0 
keV), counts versus phase, is shown in the top panel of Fig 9. 
We find some evidence for two pulses per period. Fitting the 
lightcurve with two Lorentzians (as Kuiper et al. 2002a) we 
find that the FWHM of the pulses is (50i=O.lO±O.O7, centred 
at phase 01=0.36+0.04 and ^02=0.09+0.03, centred at phase 
02=0.80+0.03. Fitting with two Gaussians gives similar results. 
Using a Z] test (Buccheri et al. 1983) on the X-ray data, we de- 
termine a value of 13, which corresponds to a probabiUty that 
the pulse-phase distribution deviates from a statistically flat dis- 
tribution of 0.99. We find a pulsed fraction of 77±13%. 

In the lower panel in Fig 9 we have plotted the ra- 
dio lightcurve from the timing observations taken using the 
Effelsberg and Lovell telescopes. 

5. Discussion 

We have investigated the two faint miUisecond pulsars PSR 
J0751 + 1807 and PSR J1012+5307 in the X-ray band 0.2-10.0 
keV, to ascertain the nature of the X-ray spectra. There are only 
a handful of MSPs that have been seen to pulse in both the ra- 
dio and in X-rays to date (Becker & Aschenbach 2002), thus 
we have also investigated whether these two MSPs show pul- 
sations in X-rays. 



PSRJ0751+1807 andPSRJ1012+5307 

Table 4. Ephemeris of PSR J1012+5307 from the Effelsberg and the 
Lovell radio timing data. Errors on the last digits are shown in paren- 
thesis after the values. 



Parameter 


Value 


Right Ascension (J2000) 


10" 12"" 33^43463677 


DecUnation (J2000) 


53" 07' 02" 4965199 


Period (P) 


0.00525605523 s 


Period derivative (P) 


0.726912 x 10-20 ^ ^-i 


Second period derivative (P) 


2.48928 x 10-30 s s^^ 


Frequency (y) 


190.267837551251880(508) Hz 


Frequency derivative (v) 


-6.200023(202)xl0-"* Hz s"' 


Second frequency derivative (v) 


1.871(294) X 10-27 Hz ^-2 


Epoch of the period (MJD) 


52018.324635450415 


Orbital period 


52243.7224464(17) s 


a.sin i 


0.581817416(121) 


Eccentricity 


<1.3 xlO-*^ 


Time of ascending node (MJD) 


52018.268144542(23) 



3.0 



2,5 - 




0.0 I 

-Z.QxlQ-* -LOxlQ-* LOxlQ-* 2.0x10-^ 

A 1/ 

Fig. 8. xl versus change in frequency from the expected pulsation 
frequency (shown as the solid vertical line at A f = 0.0) for PSR 
J1012+5307. 



Saito et al. (1997) showed that the MSP PSR B1821-24 
has a similar magnetic field value, at the light cylinder radius 
(Bl), as the Crab pulsar. If this value can be taken to indi- 
cate the high-energy magnetospheric activity, they state that 
one can expect to see similar pulses from pulsars with a similar 
value of Bl and thus such pulsars are good candidates amongst 
the MSPs in which to look for magnetospheric X-ray pulsa- 
tion in high energy bands. We have calculated the Bl for PSR 
J0751+1807 and PSR J1012+5307, using the values for the 
magnetic field strength from Zhang & Harding (2000). We 
find 8.3x10"* and 3.5x10"* G respectively. The value calculated 
for PSR J075 1 + 1 807 indicates fliat it is amongst the top 1 val- 
ues for a millisecond pulsar. That of PSR J1012+5307 places 
it in the top half. Thus if a high Bl indicates high-energy mag- 
netospheric activity this adds support to the detection of X-ray 
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Fig. 9. Upper panel: Lightcurve folded on the radio ephemeris and 
binned into 12 bins, each of 0.44 msecs. Two cycles are shown for 
clarity. A typical ±lcr error bar is shown. The dashed line shows the 
background level, where the error bar represents the ±lcr error. Lower 
panel: Radio profile from the Effelsberg and the Lovell radio timing 
data. Again two cycles are shown for clarity. 



pulsations in our observation of PSR J075 1 + 1 807 and possibly 
tothatofPSRJ1012+5307. 

We have also calculated the spin down energy (E) and the 
luminosity in the 0.1-2.4 keV band of these two pulsars, to 
compare the results with the correlation found between these 
two parameters by Becker & Triimper (1997). Becker & 
Triimper (1997) suggest that pulsars that obey this relationship 
emit X-rays produced by magnetospheric emission, originat- 
ing from the co-rotating magnetosphere. For PSR J075 1-1-1807 
we find log(E) of 34.35. To calculate the luminosity, we have 
used the distance calculated using the radio dispersion mea- 
sure and the model of Taylor & Cordes (1993). This gives a 
value of log(L;c) of 31.0+0.2, which places the point close to 
the expected value of 3 1 .2 using the correlation proposed by 
Becker & Triimper (1997). However, using the more recent 
model of Cordes et al. (2003a,b) diminishes the distance by 
almost a factor 2, to 1.1 kpc. This gives a log(L^) of 30.5+0.2 
and thus displaces the point further from the expected value. 
For the MSP PSR J1012+5307 we find a log(E) of 34.2 and 
a log(Lj:) of 30.4+0.3, where the expected log(Lj;) is approx- 
imately 31.1. This places the point quite some distance from 
the expected value, which may indicate that the relationship is 
not a hard and fast rule. Indeed Pavlov et al. (2002) found re- 



cently that when analysing pulsars in the 2.0-10.0 keV band, 
this relationship does not hold true. 

5.1. PSRJ0751+1807 

We find that the best fitting model to the X-ray spectrum of 
PSR J075 1 + 1807 is a power law. This is indicative of a mag- 
netospheric origin of the X-ray emission. We have shown ev- 
idence that there appears to be a single broad pulse emitted 
from this MSP, where the pulsation appears to change only 
slightly with increasing energy (becomes slightly broader). The 
observations also indicate that the pulsed fraction possibly in- 
creases at higher energies. A broad pulsation can be observed 
from pulsars showing hard magnetospheric emission or soft 
thermal emission and thus we can not yet discern the origin 
of the X-ray emission from PSR J0751+1807 with this data 
set. Indeed this pulsar may be best fitted by a multi-component 
model, as is the case with several other brighter MSPs e.g. PSR 
J0218+4232, (Webb et al. 2003) or PSR J0437-4715, (Zavlin 
& Pavlov 1998; Zavlin et al. 2002). Longer observations of 
PSR J075 1 + 1 807 will help to distinguish the true nature of the 
spectrum. 

5.2. PSR J1 01 2+5307 

We can not discriminate which of the model fits to the MSP 
PSR J1012+5307 is the best, however we find that the single 
power law has a similar photon index to the X-ray spectrum of 
PSR J0751+1807 (see Sect. 3.1), which could indicate a mag- 
netospheric origin of the X-ray emission. However, the tem- 
perature of the blackbody (1.9+0.5 xlO^ K) is consistent with 
that emitted from the heated polar caps of a millisecond pul- 
sar (10^-10^ K e.g. Zhang & Cheng 2003; Zavlin & Pavlov 
1998, and references therein). Calculating the radius of the 
emission area from the blackbody model fit, we find a radius of 
0.05+" " km, which is smaller than the expected radius of emis- 
sion from polar caps (~1 km e.g. Zhang & Cheng 2003; ZavUn 
& Pavlov 1998, and references therein). Zavlin et al. (1996) 
state however, that spectral fits with simplified blackbody mod- 
els can produce higher temperatures and smaller sizes due to 
the fact that the X-ray spectra emerging from light-element 
atmospheres are harder than blackbody spectra. Alternatively, 
Zavlin & Pavlov (1998) and Zavlin et al. (2002) suggest that 
the thermal emission can be from non uniform polar caps and 
we may therefore be seeing the emission from the hotter central 
region of the caps. 

Folding the timing data on the radio frequency, we find 
some evidence for two pulses emitted from this MSP, separated 
by approximately 0.4 in phase. This is similar to other millisec- 
ond pulsars, e.g. PSR J0218+4232 (Webb et al. 2003) and PSR 
B1821-24(Saitoetal. 1997). 

6. Conclusions 

We have presented XMM-Newton data of the faint millisec- 
ond pulsars PSR J0751+1807 and PSR J1012+5307. Both of 
these pulsars have a reasonably large magnetic field at the light 
cylinder radius, which could indicate that both of these MSPs 
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should show pulsations in X-rays. We present for the first time 
the X-ray spectra of these two faint millisecond pulsars. We 
find that a power law model best fits the spectrum of PSR 
J0751 + 1807, r=1.59±0.20, with an unabsorbed flux of 4.4 
xlO'^'^ergs cm"^s"^ (0.2-lO.OkeV). Apowerlawisalsoagood 
description of the spectrum of PSR J1012+5307, r= 1 .78+0.36, 
with an unabsorbed flux of 1.2 xlO '-^ergs cm^^s ' (0.2-10.0 
keV). However, a blackbody model can not be excluded as 
the best fit to the data. We have also shown evidence to sug- 
gest that both of these MSPs may show X-ray pulsations. PSR 
J075 1+1807 appears to show a single pulse, whereas PSR 
J1012+5307 may show some evidence for two pulses per pulse 
period. 
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